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Abstract: We investigate the nonlinear response of photonic crystal 
waveguides with suppressed two-photon absorption. A moderate decrease 
of the group velocity (~ c/6 to c/15, a factor of 2.5) results in a dramatic (x 
30) enhancement of three-photon absorption well beyond the expected scal- 
ing, oc 

l/v\. This non-trivial scaling of the effective nonlinear coefficients 
results from pulse compression, which further enhances the optical field 
beyond that of purely slow-group velocity interactions. These observations 
are enabled in mm-long slow-light photonic crystal waveguides owing to 
the strong anomalous group- velocity dispersion and positive chirp. Our 
numerical physical model matches measurements remarkably. 
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1. Introduction 

Slow-light nonlinearities have been remarkably observed through quantum coherence and inter- 
ference in atomic systems, with group velocities of tens of meters per second and sub- 100 kHz 
bandwidths d [2 EJ |4j . In solid-state systems such as photonic crystal cavities or waveguides 
modest slow group velocities of down to c/300 in comparison has been achieved but 
possess THz bandwidths for chip-scale optical signal processing OHO). The strong ab-initio 
structural dispersion in slow-light photonic crystals waveguides (PhCWGs) not only gives rise 
to localized modes HU [121 Q2L but also dramatic enhancement of resonant and non-resonant 
nonlinearities lfT4l [T5l . Third-order processes have also recently been observed for slow-light 
fT6lL self-phase modulation (SPM) El ED [19), anc * third-harmonic generation in optical mi- 
crostructures O EH [22) . Here we present the first observations of self-phase modulation 
limited only by three-photon- absorption as well as evidence for pulse compression in slow- 
light GalnP photonic crystal waveguides. In contrast to previous self-phase modulation work 
in PhCWGlT3l [T8l Q2), which was limited by the third-order two-photon absorption (TPA), 
here the nonlinear loss term is the fifth-order three-photon absorption. The different orders of 
the desirable and undesirable nonlinearities, along with tight modal confinement of the pho- 
tonic crystal, lead to the observation of novel optical effects. Remarkably the three-photon 
absorption process demonstrates a 30-fold enhancement, and departs from the expected (1/v^) 
scaling, even when taking into account slow-light disorder scattering. Pulse compression was 
further observed in our positive chirp due to the interaction of the Kerr effect with anomalous 
group velocity dispersion(GVD). These series of measurements show excellent match with our 
numerical simulations including group-velocity-dependent nonlinearities and losses, and meas- 
ured group velocity dispersion through optical low-coherence reflectometry l23l . 

The origin of slow-light in photonic crystals arises from coherent Bragg reflections due to 
the in-plane periodic PhCWG lattice, leading to an ultra-flat dispersion of the transverse electric 
(TE) field. At these so-called slow-light frequencies, the light effectively travels slowly through 
the lattice via multiple Bragg reflections, leading to an enhanced local field density. The local 
field enhancement scales inversely with the group velocity, thus decreasing the threshold of 
intensity-dependent nonlinear effects such as Kerr, multi-photon absorption, or Raman scatte- 
ring O [141 Q2) . The combination of the engineered group velocity enhancement, along with 



the small modal effective area of PhCWGs (^2x 10 -13 m 2 ), yields outstanding control over 
the optical modes. 

The optical Kerr-effect has long been studied for its utility in nonlinear all-optical de- 
vices (l4l[24]|. In addition to the desirable Kerr term, each material also has a nonlinear multi- 
photon absorption term. This fundamental material property induces several detrimental effects, 
such as: (a) limiting the Kerr-induced phase- shift, also called self-phase modulation (SPM); (b) 
inducing nonlinear and free-carrier absorption losses; (c) distorting the pulse via free-carrier 
dispersion; as well as (d) restricting the spectral range of any potential nonlinear optical devices 
such as all-optical switches l25l . In fact, several early studies point out that the wavelength range 
of the two-photon absorption (TPA) tail depends strongly on the quality of the molecular beam 
epitaxy (MBE) sample growth [25, 26,23- These fundamental obstacles must be overcome in 
order to achieve practical nonlinear devices, such as all-optical switches, on-chip. With this in- 
sight, we selected the material system GalnP (E g =l.9 eV) for its desirable nonlinear properties. 
With an energy bandgap at least 300 meV above the TPA range (E g > 2hco), the GalnP sample 
investigated here is well outside the range of potential defect states and the sole nonlinear loss 
mechanism is three-photon absorption (ThPA). Importantly, ThPA is negligible at the optical 
intensities required for the Kerr-effects (> n phase-shifts) observed experimentally in this pa- 
per. The complete suppression of TPA and small impact of ThPA on the GalnP PhCWG open 
the possibility of a significant spectral window for all-optical signal processing (251 . 
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Fig. 1. (Color online) Linear properties: (a) Measured PhCWG transmission and corre- 
sponding n g of 6.6, 7.5, 8.8, 10.7, 12.6, 14.2, respectively, from low-coherence reflectome- 
try after Refs. I231I3T1 . Inset: SEM image with scale bar of 1 jUm. (b) Sample optical low- 
coherence reflectometry (OLCR) data used to extract n g = c^j^. Inset: Extracted group 
indices versus derived GVD coefficients. 



2. Experiment 



2.1. Sample details and linear characterization 

Our PhCWG sample is a Wl GalnP membrane of 190 nm thickness, with a hexagonal lattice 
constant a of 480 nm, hole radius r w 0.23a, and 1.5 mm length. The fabrication has been 
described elsewhere [[28l and the linear PhCWG transmission is illustrated in Figure [Ha). We 
carefully designed integrated mode-adapters |[29) to reduce losses to ~ 5 dB/facet at 1526 
nm BUI and suppress Fabry-Perot oscillations, though disordered propagation due to the lattice 
is clearly seen at the slow-light onset region ll3Tll . Optical low-coherence reflectometry [[2311371 
[Fig- Eh)] was used to extract the group indices, and the group velocity dispersion (GVD) 
[Fig- Eh) inset] computed using complete 3D planewave expansion |[32lL with measurement- 
consistent group indices. 

2.2. Nonlinear characterization - three-photon absorption 

For the nonlinear measurements, we employed a wavelength-tunable mode-locked fiber laser 
(PriTel) with w 3 to 5 ps pulses [[33l (characterized with autocorrelator) at 22 MHz with 
transverse-electric (electric-field in-plane) polarization. The output pulse monitored with a 
spectrum analyzer and an oscilloscope. The peak power coupled to the PhCWG P c is defined 
as: P c = T]K c Pi n , where the input peak power P[ n (output power P out ) includes the objective loss 
(2 dB), mode mismatch (2 dB) r\ |29lL with the coupling coefficient, K c . The insertion loss due 
to disorder varies with slow-light 1341 : the values of which are extracted directly from Eqn.Q] 
below. 

In Fig. [2a), we first examined P[ n — P out at different group velocities and same input power 
range to illustrate the enhanced nonlinear absorption processes. At the largest group velocities, 
Pin — Pout is linear while at smaller group velocities the output power begins to saturate, indicat- 
ing the distinct onset and slow-light enhancement of ThPA, an intensity dependent loss mech- 
anism, in the slow-light regime. The nonlinear propagation equation (dP/dz = — ocP — CC3P 3 ) 
with symmetric input-output coupling gives: 

_L = JjL = 1 , 2 ^effL3eff 2 

T 2 Plut n 2 K*e- 2aL ^ K 2 e~ 2 ^ in " K) 

where OL^ e ff — jp - {jt) includes the 5th-order modal effective area, and L^ e ff = (1 — 

e~ 2aL )/(2a). The bulk 0C3 coefficient, ~6x 10 _26 m 3 /W 2 , was calculated from a well-cited 
model[ 35], experimentally verified for the similar AlGaAs material in Ref. ll27l . 

In Fig. Ob), we illustrate an example plot of 1 /T 2 versus Pf n to extract the effective nonlinear 
ThPA coefficient and coupling factor K c at a particular group index n g of 8.8. We note the 
material-dispersion of the nonlinear susceptibilities is negligible within our measurement 
range. An example two-photon absorption (l/T versus Pi n ) analysis l27l [FigEJc)] shows 
clear mismatch of the experimental data (blue squares) with the TPA fit (dashed line). Similar 
fits occur for all group velocities investigated, thus confirming the sample experiences only 
three-photon absorption. We note that, for the group velocities we examined, the linear loss 
coefficient a has been observed to scale |[T2l l34l l36l approximately as n 2 , with our measured 
values of 1 dB/mm at 1526 nm BTIl . In the semilog plot of Fig. [2d), we illustrate the extracted 
three-photon (Xn> e ff values in the slow-light regime. In this first observation of three-photon 
absorption in PhCWGs, we further emphasize that a surprising 30-fold enhancement in 
three-photon coefficient from 0.15 to 4 W~ 2 cm~ l was observed with n g increasing only from 
6.6 to 14.2. While the field enhancement of nonlinear absorption due to slow-light can be 
predicted from the material^) and device parameters (n g ) and is known to obey a prescribed 




Fig. 2. (Color online) Nonlinear measurements: (a): P\ n —P ou t depending on the group ve- 
locity. The increased curvature corresponds to larger nonlinear absorption due to slow-light 
enhanced three-photon absorption (ThPA) at longer wavelengths. Key: 1526 nm(circles), 
1530nm(diamonds), 1534 nm(squares), 1538 nm(upward triangles), 1541 nm(downward 
triangles), and 1544 nm(stars). (b) Sample plot of the inverse transmission squared (l/T 2 ) 
versus Pf n at 1534 nm (n g =S.S). The points are experimental data and the line is the best 
fit to extract the effective ThPA coefficient (%3eff- The key is the same as in (a), (c) Exam- 
ple inverse transmission (l/T) versus P c plot depicting the mismatch of (and negligible) 
two-photon absorption in our slow-light GalnP PhCWGs, for the same experimental data 
of panel (b). (d) Extracted CL?,eff (black dots) versus group index with the expected scaling 
of ThPA (solid red curve) 



trend lfT4l . here the deviation in the shape of the trajectory with group index (n g ) is markedly 
different than the scaling predicted from exclusively slow-light effects, an indication that 
some other field enhancement effect is involved. Moreover, this enhancement is by far faster 
than the predicted scaling from the definition of the ThPA coefficient l38l , — ) 3 ~ 5.8 

A 5eff U ° 

enhancement expected from simple slow-light effects. We now further elaborate on this point 
in terms of pulse self-phase modulation due to the optical Kerr effect. 



2.3. Nonlinear characterization - Kerr-effect 

To further illuminate this nontrivial scaling, we measured the output spectra at different slow 
group velocities [1526, 1534, and 1538 nm in Figs.[3ta), (b), and (c) respectively ]. At increased 
input intensities, the pulses undergo stronger Kerr self-phase modulation with the associated 
spectral broadening. In each plot, we show examples of low (dot-dashed black line) and high 
(>2 W; solid blue line) power spectra. For group indices less than ~ 10, a clear increase in 
spectral width was observed with increasing group index, while maintaining the same coupled 
peak power. The fine-peaked structure in the pulse spectra is due to disorder-enhanced scatte- 
ring |3H . The high-power spectra exhibit symmetric double-peak structures characteristic of 



Kerr-induced SPM, where the pulse spectral symmetry rules out free-carrier dispersion gener- 
ated by ThPA |QJ) or third-order dispersion B91I401 . 



Table 1 . Parameters used in numerical simulations 



Parameter 


Value 


Value 


Value 


Value. 


Wavelength X (nm) 


1526 


1530 


1534 


1538 


Pulse duration (FWHM) r(ps) 


5.0 


4.4 


4.5 


3.0 


Chirp C 


1.3 


1.1 


1.2 


0.2 


Time-bandwidth product Av T 


0.74 


0.69 


0.69 


0.33 


Peak Power P C (W) 


1.85 


2.2 


2.25 


2.35 


Group index n g 


6.6 


7.5 


8.8 


10.7 


GVD p 2 (ps 2 /mm) 


-0.6 


-0.97 


-1.4 


-2.2 


Eff. ThPA a 3ieff (cm- l W- 2 ) 


0.24 


0.40 


0.61 


0.66 


SPM gain yicm^W' 1 ) 


6.75 


7.96 


9.7 


13.0 



To further quantify our experimental results, we performed numerical simulations of the 
nonlinear Schrodinger equation |4H : 

where A is the pulse envelope amplitude with P c = |A| 2 , ft [ps 2 /mm] the GVD parameter, z is 
the propagation direction, and t defined in a moving reference frame. The effective nonlinear co- 
efficients are defined according to Ref. (42l|43) and modes from planewave expansion l32l . We 
computed A$ e ff for the fifth-order nonlinearity, ThPA, and A 3e ff for the third-order Kerr nonlin- 
earity. Our calculations show that A$ e ff & 0.15A 3e ff over the wavelength range and waveguide 
geometry examined here. The third-order area A^eff ranged from 2 and 2.6 x 10 -13 m 2 for the 
wavelengths examined, an increase far smaller than the slow-light scaling. Assuming a conser- 
vative absorbed energy per pulse due to ThPA of 20%, for a (~4 ps) this corresponds to ~pJ 
or less distributed along the full length of the waveguide (1.5 mm). At the repetition rate (22 
MHz) of the laser, the average absorbed power would then be on the order of tens of /iW. Thus 
at the power levels used in the experiment, free-carrier effects contribute negligibly to the pulse 
dynamics in the GalnP material. Positive chirped (chirp parameter C) hyperbolic- secant input 
pulses |'4l], such as those generated by the mode-locked fiber laser in the experiment, are used. 
The variation of the chirp over the wavelength range is a characteristic property of the laser. 
A Kerr n 2 coefficient of 8 x 10 _14 cm 2 /W l44lh scaled with n 2 , along with the experimentally- 
extracted values of OL^eff are use d in our model. Table Q] summarizes the measured parameters 
in our model, with the peak power within 10% of measured values (uncertainty from chirp and 
pulsewidth). 

The resulting simulated spectra with chirp and GVD at three different group indices, ft g =6.6, 
8.8, and 10.7 in Figs. [2a), (b), and (c), respectively, show remarkable matches with our 
measurements. To examine the effects of chirp and GVD on ThPA, in Fig. Od), (e), and (f) 
we demonstrate 1/T 2 versus P 2 generated strictly from experimental parameters. The complete 
physical model not only describes the higher-order nonlinear absorption (1/T 2 versus P 2 ), but 
also captures its scaling to larger values in both slope and curvature at lower group velocities. 
The inclusion of anomalous GVD in the model, and thus the possibility of pulse compression, 
not only reproduces the slight upward bending, but also rigorously and correctly predicts the 
increase of the effective absorption beyond the n 3 g scaling. We note that higher order dispersion 
is still negligible. The slight upward bending in 1/r 2 vs. P 2 n plot is the signature of pulse- 
compression. We emphasize that despite this small bending, the fit for (X^eff are quite good for 




Fig. 3. (Color online) Comparison of experimental and theoretical spectra for different 
wavelengths of 1526 nm (n g =6.6) (a), 1534 nm (n g =S.S) (b), and 1538 nm 0^=10.7) (c), 
with chirped sech 2 input. Experimentally derived parameters, also used in simulations, are 
shown in Table [T] (d) Experimentally l/T 2 versus P 2 for 1526 nm(n g =6.6). The solid 
line shows the simulated results including GVD, while the dashed line are simulation re- 
sults without GVD. (e),(f) Same as (c) for but with (e) 1534 nm(^=8.8) and (f) 1538 
nm(n g =\0.1). The strong upward bend of the curve indicates enhanced nonlinear absorp- 
tion beyond conventional slow-light scaling, triggered from pulse compression. 



all values examined. We also examined the case where GVD is set to zero as the dotted line in 
FigS-Etd-f). The strong deviation between the simulation when GVD is zero and experimental 
data clearly demonstrate that GVD is interacting strongly with the nonlinear effects. We also 
simulated the absence slow-light effects, that is n g = 3.12. In all cases, the curves were well 
below the experimental data, and nearly overlapped the zero GVD simulation data. These are 
not included here for figure clarity. 

To further examine the impact of pulse compression on the slow-light Kerr nonlinearity, we 

quantified the measured SPM-induced spectra broadening with: AA 2 = ^j^p^ ' w h ere S 
is the lineshape. Illustrative measurements of spectral broadening at n g =6.6, 8.8, and 10.7 for 
increasing coupled power are shown in Fig. Ufa), (b), and (c), respectively. While SPM alone 
(dashed line) produces a linear slope that has a slight downward curvature due to higher-order 
absorption for large n g , the spectral widths including the effects of GVD mark a significant 
departure from a linear slope, as observed both experimentally (markers) and numerically(solid 
line). Fig.^c) in particular demonstrates a saturation of the spectral broadening at larger in- 
put powers. This is confirmed to arise from a combination of ThPA and GVD in the numeri- 
cal model. We additionally plot the case without slow-light effects as the dashed line. Clearly 
slow-light has a strong impact on the spectral broadening, as the behavior is completely linear. 
Similar results occur at the other wavelengths. 
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Fig. 4. (Color online)(a) and (b) RMS pulse broadening (nm) as a function of coupled 
input power P c at: (a) 1526 nm (n g =6.6), (b) 1534 nm (n g =S.S) and (c) 1538 nm (n g =\0J). 
The points are experimental data, the solid line is simulation with GVD, and the dashed 
line is without GVD. We also show the case without slow-light (e.g. n g = 3.12) in (c) as 
the dashed line, (d) N = ^Ld~/Lnl versus n g at P[ n — max(W), the max peak power input 
into the PhCWG. For values of N > 1, the pulse has the possibility of being compressed, 
(e) Plot of L/Lp vs. rig. In addition to N > 1, the pulse must also propagate a minimum 
length, related to the dispersion length L#, before compression can occur, (f) Effective 
nonlinear absorption, CL^eff an d square of the effective SPM coefficient, y^y, rescaled with 
a suitable constant, C. The local field enhancement of the two effects scales as predicted. 
The experimental values demonstrate non-trivial scaling. 



2.4. Discussion of non-trivial scaling 

Thus far we have observed trends in both our ThPA and Kerr measurements that do not 
obey conventional slow-light scalings. Clearly a field enhancement effect beyond slow-light 
enhancement is occurring in the sample. Based on the necessity of including GVD in our 
modeling, and with knowledge that pulse compression (resulting from the interaction of the 
Kerr self-phase modulation and anomalous dispersion) causes increased peak power of an 
optical pulse flTTl . we investigated pulse compression as a source of the additional local field 
in our PhCWG. We first computed the dispersion length, Ld, (~ 14 mm at 1526 nm and 
1.3 mm at 1538 nm) and the nonlinear length, Lnl, (~ 1.3 mm at 1526 nm and 0.45 mm 
at 1538 nm) BTll at the various pulse center frequencies of our slow-light experiments. The 
results, N = ^Ld/Lnl, are plotted in Fig. HJd). For N > 1 and ftC < 0, pulse compression 
is possible[41]. There is, however, an additional criterium to observe pulse compression. 
The requirement is that the pulse must propagate a minimum distance compared to Lp, the 
dispersion length, and depends on N. This implies that the physical device length must be 
comparable to Ld, typically at least Ld/2 BTll . Figure HJe) demonstrates the second key 
parameter, the L/Ld ratio, for the experimental data. While N is always greater than one, the 
dispersion length here decreases such that L/Ld increases from L/Ld = 0.1 to Ld = 1.4 as n g 
goes from 6.6 to 12.6. Thus while the pulses of the lowest three group indices n g experience 



some compression (N > 1 and L/Lj) « 0.1 — 0.3), the larger ftg values experience more 
compression since they possess both N > 1 and ratios of L/L^ > 0.5. Thus both criteria have 
been met. This accounts for the faster than predicted scaling of Cfo>eff- 

We now utilize the extracted broadening data to determine the effective SPM coefficient Yeff 
from the nonlinear phase (j) max = YeffPc^eff and spectral broadening: ^ = + ^/^Qmax- 
Here the effective length L e ff is limited only by ThPA and hence must be defined differently 
from the typical linear loss parameter l45l . The enhancement of the three-photon absorption 
and SPM broadening are related through the field enhancement factor, k, relating the propa- 
gating pulse power to the electric field energy density: \E\ 2 = kP. We thus have Yeff ^ K an d 
a 3,eff ^ 2 , implying that Oh> e ff scales as Yeff- In Fig.^f), we examine Oh> e ff (black circles) 
and y^ff ( re d squares; scaled by a constant C) - the curves are remarkably similar, highlighting 
unambiguously that the non-trivial scalings of Kerr and ThPA have a common physical origin, 
that of strong field enhancement consistent with pulse compression. 

3. Conclusion 

We have demonstrated nonlinear scaling of the nonlinear enhancement beyond the limits 
of slow-light photonic crystal waveguides. This non-trivial scaling has been analyzed with 
experimental data of both self-phase modulation and three-photon absorption and further 
reinforced with a nonlinear propagation model. The origin of this non-trivial scaling beyond 
the slow-light regime is derived from an additional local field enhancement due to pulse 
compression. As Kerr is a third-order effect in field, while ThPA is a weaker fifth-order 
nonlinear effect, dispersion engineering to include pulse compression can greatly reduce the 
threshold of SPM, while carefully managing ThPA and its undesirable properties. The ability 
to engineer PhCWG dispersion! 8 1, such as with low group velocity dispersion [|46l to precisely 
control the pulsewidth, while completely suppressing TPA, make GalnP PhCWGs promising 
for compact broadband ultra-fast optics. 
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